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Abstract 

We extract the e + e _ —> 7r + 7r - cross section in the energy range between 600 and 900 MeV, exploiting the 
method of initial state radiation. A data set with an integrated luminosity of 2.93 fb _1 taken at a center- 
of-mass energy of 3.773 GeV with the BESIII detector at the BEPCII collider is used. The cross section 
is measured with a systematic uncertainty of 0.9%. We extract the pion form factor |F„-| 2 as well as the 
contribution of the measured cross section to the leading-order hadronic vacuum polarization contribution 
to (g — 2)^. We find this value to be a^ 7r,LO (600 — 900 MeV) = (368.2 ± 2.5 sta t ± 3.3 sys ) • 10~ 10 , which is 
between the corresponding values using the BaBar or KLOE data. 

Keywords: Hadronic cross section, muon anomaly, initial state radiation, pion form factor, BESIII 


1. Introduction 

The cross section aw = a(e + e~ —> 7r + 7r~) 
has been measured in the past with ever increas¬ 
ing precision at accelerators in Novosibirsk aim 
Orsay [4], Frascati !H[6l HUE], and SLAC [H Ho]. 
More recently, the two most precise measurements 


have been performed by the KLOE collaboration in 
Frascati [Bj and the BaBar collaboration at SLAC 
mm- Both experiments claim a precision of bet¬ 
ter than 1% in the energy range below 1 GeV, in 
which the p(770) resonance with its decay into pions 
dominates the total hadronic cross section. A dis¬ 
crepancy of approximately 3% on the peak of the 
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p( 770) resonance is observed between the KLOE 
and BaBar spectra. The discrepancy is even in¬ 
creasing towards higher energies above the peak of 
the p resonance. Unfortunately, this discrepancy 
is limiting the current knowledge of the anomalous 
magnetic moment of the muon = (<7 — 2)^/2 [TT] . 
a precision observable of the Standard Model (SM). 
The accuracy of the SM prediction of (g — 2) M is en¬ 
tirely limited by the knowledge of the hadronic vac¬ 
uum polarization contribution, which is obtained in 
a dispersive framework by using experimental data 
on a(e + e~ — > hadrons) [Til 1T2I fT5] , The cross sec¬ 
tion a(e + e~ —> 7 r + 7 r _ ) contributes to more than 
70% to this dispersion relation and, hence, is the 
most important exclusive hadronic channel of the 
total hadronic cross section. Currently, a discrep¬ 
ancy of 3.6 standard deviations (T2J is found be¬ 
tween the direct measurement of and its SM pre¬ 
diction. However, the discrepancy reduces to 2.4cr 
0 . when only BaBar data is used as input to the 
dispersion relation. In this letter we present a new 
measurement of the cross section ov^, obtained by 
the BESIII experiment at the BEPCII collider in 
Beijing. 

The measurement exploits the method of initial 
state radiation (ISR), the same method as used by 
BaBar and KLOE. In the ISR method events are 
used in which one of the beam particles radiates 
a high-energy photon. In such a way, the avail¬ 
able energy to produce a hadronic (or leptonic) fi¬ 
nal state is reduced, and the hadronic (or leptonic) 
mass range below the center-of-mass (cms) energy 
of the e + e~ collider becomes available. In this pa¬ 
per, we restrict the studies to the mass range be¬ 
tween 600 and 900 MeV/c 2 , which corresponds to 
the p peak region. 

The remainder of this letter is organized as fol¬ 
lows: In section 2, the BESIII experiment is in¬ 
troduced. In section 3 we describe the data set 
used, the Monte Carlo (MC) simulation, the event 
selection of e + e _ —> 7 r + 7 r - 7 events, and the data- 
MC efficiency corrections. The determination of the 
integrated luminosity of the data set is described 
in Section [4] A cross check of the used efficiency 
corrections using the well-known e + e _ —> p + p~"/ 
QED process is performed in Section [5j before ex¬ 
tracting the 7 r + 7 r^ cross section in Section [ 6 j 

2. The BESIII experiment 

The BESIII detector is located at the 
double-ring Beijing electron-positron collider 


(BEPCII) ng. 

The cylindrical BESIII detector covers 93% of 
the full solid angle. It consists of the following de¬ 
tector systems. (1) A Multilayer Drift Chamber 
(MDC), filled with helium gas, composed of 43 lay¬ 
ers, which provides a spatial resolution of 135 pm, 
an ionization energy loss dE/dx resolution better 
than 6 %, and a momentum resolution of 0.5% for 
charged tracks at 1 GeV/c. (2) A Time-of-Flight 
system (TOF), built with 176 plastic scintillator 
counters in the barrel part, and 96 counters in the 
endcaps. The time resolution is 80 ps in the barrel 
and 110 ps in the endcaps. For momenta up to 1 
GeV/c, this provides a 2 a K/ 7 T separation. (3) A 
CsI(Tl) Electro-Magnetic Calorimeter (EMC), with 
an energy resolution of 2.5% in the barrel and 5% 
in the endcaps at an energy of 1 GeV. (4) A su¬ 
perconducting magnet producing a magnetic field 
of IT. (5) A Muon Chamber (MUC) consisting of 
nine barrel and eight endcap resistive plate chamber 
layers with a 2 cm position resolution. 

3. Data sample, event selection, and effi¬ 
ciency corrections 

3.1. Data sample and MC simulations 

We analyze 2.93 fb _1 (see Sect. [4j) of data 
taken at a cms energy y/s = 3.773 GeV, which 
were collected in two separate runs in 2010 and 
2011. The Phokhara event generator ns nzi is 
used to simulate the signal process e + e _ —>■ 7 r + 7 r _ 7 
and the dominant background channel p + p~ 7 . 
The generator includes ISR and final state radi¬ 
ation (FSR) corrections up to next-to-leading or¬ 
der (NLO). Effects of ISR-FSR interference are in¬ 
cluded as well. The continuum qq (q = u, d, s) MC 
sample is produced with the KKMC event genera¬ 
tor [18]. Bhabha scattering events are simulated 
with BABAYAGA 3.5 [T9]. The Bhabha process is 
also used for the luminosity measurement. All MC 
generators have been interfaced with the Geant4- 
based detector simulation [ 20 ] [ 21 ]. 

3.2. Event selection 

Events of the type e + e^ —> 7 r + 7 r - 7 are selected. 
Only a tagged ISR analysis is possible in the 
mass range 600 < m V7r < 900 MeV/c 2 , where m^ 
is the 7 r + 7 r - invariant mass, i.e., the radiated pho¬ 
ton has to be explicitly detected in the detector. 
For untagged events, the photon escapes detection 
along the beam pipe; the hadronic system recoiling 


4 


against the ISR photon is therefore also strongly 
boosted towards small polar angles, resulting in 
no geometrical acceptance in the investigated m W7r 
range. 

We require the presence of two charged tracks 
in the MDC with net charge zero. The points of 
closest approach to the interaction point (IP) of 
both tracks have to be within a cylinder with 1 cm 
radius in the transverse direction and ±10 cm of 
length along the beam axis. For three-track events, 
we choose the combination with net charge zero for 
which the tracks are closest to the IP. The polar an¬ 
gle 9 of the tracks is required to be found in the fidu¬ 
cial volume of the MDC, 0.4rad < 9 < n — 0.4 rad, 
where 9 is the polar angle of the track with respect 
to the beam axis. We require the transverse mo¬ 
mentum p t to be above 300 MeV/c for each track. 
In addition, we require the presence of at least one 
neutral cluster in the EMC without associated hits 
in the MDC. We require a deposited energy above 
400 MeV. This cluster is then treated as the ISR 
photon candidate. 

The radiative Bhabha process e + e~ —» 
e + e _ 7 ( 7 ) has a cross section which is up to three 
orders of magnitude larger than the signal cross sec¬ 
tion. Electron tracks, therefore, need to be sup¬ 
pressed. An electron particle identification (PID) 
algorithm is used for this purpose, exploiting infor¬ 
mation from the MDC, TOF and EMC [22]. The 
probabilities for being a pion P( n) and being an 
electron P(e) are calculated, and P{n) > P(e) is 
required for both charged tracks. 

Using as input the momenta of the two selected 
track candidates, the energy of the photon can¬ 
didate, as well as the four-momentum of the ini¬ 
tial e + e _ system, a four-constraint (4C) kinematic 
fit enforcing energy and momentum conservation 
is performed which tests the hypothesis e + e _ —► 
7 r + 7 r - 7 . Events are considered to match the hy¬ 
pothesis if they fulfill the requirement xlc < 60- 
It turns out that the p + p~ 7 final state can not be 
suppressed by means of kinematic fitting due to the 
limited momentum resolution of the MDC. An in¬ 
dependent separation of pion and muon tracks is 
required. 

We utilize a track-based muon-pion separation, 
which is based on the Artificial Neural Network 
(ANN) method, as provided by the TMVA pack¬ 
age [23]. The following observables are exploited 
for the separation: the Zernicke moments of the 
EMC clusters [2 23, induced by pion or muon tracks, 
the ratio of the energy E of a track deposited 


in the EMC and its momentum p measured in 
the MDC, the ionization energy loss dE/dx in the 
MDC, and the depth of a track in the MUC. The 
ANN is trained using 7 r + 7 r _ 7 and p + p ~7 MC sam¬ 
ples. We choose the implementation of a Clermont- 
Ferrand Multilayer Perceptron (CFMlp) ANN as 
the method resulting in the best background rejec¬ 
tion for a given signal efficiency. The output likeli¬ 
hood 2 M.NN is calculated after training the ANN for 
the signal pion tracks and background muon tracks. 
The response value ?/ann is required to be greater 
than 0.6 for each pion candidate in the event selec¬ 
tion, yielding a background rejection of more than 
90% and a signal loss of less than 30%. 

3.3. Efficiency corrections 

Given the accuracy of 0(1%) targeted for the 
cross section measurement, possible discrepancies 
between data and MC due to imperfections of the 
detector simulation need to be considered. We have 
investigated data and MC distributions concerning 
the tracking performance, the energy measurement, 
and the PID probabilities, both for the electron PID 
as well as the pion-muon separation. In order to 
produce test samples of muon and pion tracks over 
a wide range in momentum/energy and polar an¬ 
gle, we select samples of p + p~j and 7 r + 7 r _ 7 r + 7 r _ 7 
events that have impurities at the per mille level. 
By comparing the efficiencies found in data with 
the corresponding results found in the MC sam¬ 
ples, we determine possible discrepancies. Corre¬ 
sponding correction factors are computed in bins of 
the track momentum or energy and the track po¬ 
lar angle 9 , and are applied to MC tracks to ad¬ 
just the reconstructed number of events. While 
for the reconstruction of charged tracks and neu¬ 
tral clusters and for electron PID, the differences 
between data and MC are smaller than 1 % on aver¬ 
age, differences up to 10% occur in the ANN case. 
The corrections are applied separately for neutral 
clusters and for muon and pion tracks. Hence, we 
do not only obtain the corrections for the 7 r + 7 r ~7 
signal events, but also for the dominating p + p~ 7 
background. The statistical errors of the correction 
factors are included in the statistical uncertainty 
of the measurement. Systematic uncertainties as¬ 
sociated to the correction factors are presented in 
Sect. |6.5| The efficiency correction for the photon 
efficiency is obtained after the application of the 
kinematic fit procedure. The corresponding correc¬ 
tion is therefore a combined correction of photon 
efficiency and differences between data and MC of 
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the xlc distribution. The systematic uncertainty 
for the contribution of the photon efficiency and 
x\c distribution is, hence, incorporated in the sys¬ 
tematic effects associated with the efficiency correc¬ 
tions. The systematic uncertainty connected with 
the pt requirement is also associated with the cor¬ 
responding efficiency correction. 

3-4- Background subtraction 

The p + p~"f background remaining after the ap¬ 
plication of the ANN is still of the order of a few per¬ 
cent, compared to 5x10 s signal events. It is, how¬ 
ever, known with high accuracy, as will be shown 
in the next section, and is subtracted based on MC 
simulation. Additional background beyond p + 
remains below the one per mille level. Table 1 lists 
the remaining MC events after applying all require¬ 
ments and scaling to the luminosity of the used data 
set. 

Table 1: Total number of remaining non-muon 
background events between 600 < m< 900 
MeV/c 2 obtained with MC samples. 


Final state 

Background events 

C^ 

1 

<D 

+ 

12.0 ±3.5 

7T + 7T 7T°7 

3.3 ± 1.8 

7r + 7r _ 7r°7r°7 

negl. 

K+K-'y 

2.0 ± 1.5 

K°K^'y 

0.4 ±0.6 

PP1 

negl. 

continuum 

3.9 ± 1.9 

ip( 3770) -> D+D- 

negl. 

ip( 3770) ->• D°D° 

negl. 

ip(3770) —> non DD 

3.1 ± 1.8 

7 ip(2S) 

negl. 

7 J/ip 

0.6 ±0.8 


4. Luminosity measurement using Bhabha 

events 

The integrated luminosity of the data set used 
in this work was previously measured in Ref. 1241 
with a precision of 1.0% using Bhabha scattering 
events. In the course of this analysis, we re-measure 
the luminosity and decrease its systematic uncer¬ 
tainty by the following means: (1) Usage of the 
BABAYAGA@NLO [25] event generator with a the¬ 
oretical uncertainty of 0.1%, instead of the previ¬ 
ously used BABAYAGA 3.5 event generator with an 


uncertainty of 0.5% m- ( 2 ) Precise estimation of 
the signal selection efficiencies. In particular, the 
uncertainty estimate of the polar angle acceptance 
is evaluated by data-MC studies within the fiducial 
EMC detection volume, which is relevant for the 
luminosity study (0.13%). The very conservative 
estimate in [ 21 ] was based on acceptance compar¬ 
isons with and without using the transition region 
between the EMC barrel and endcaps, leading to 
additional data-MC differences (0.75%). The other 
uncertainties of [241 remain unchanged and addi¬ 
tional systematic uncertainties due to the uncer¬ 
tainty of yfs ( 0 . 2 %) and the vacuum polarization 
correction (< 0.01%) are taken into account. Fi¬ 
nally, the total integrated luminosity amounts to 
C = (2931.8 ± 0.2 sta t ± 13.8 sys )pb _1 with a relative 
uncertainty of 0.5%, which is consistent with the 
previous measurement [24] . 

5. QED test using e+e — —> events 

The yield of events of the channel 
e + e _ —> n + fi~'y as a function of the two-muon 
invariant mass can be compared to a pre¬ 

cise prediction by QED, which is provided by the 
Phokhara generator. We select muon events ac¬ 
cording to the ANN method described previously 
and require ?/ann < 0.4 for both tracks, resulting in 
a background rejection of more than 90% and a sig¬ 
nal loss of less than 20%. All other requirements in 
the selection are exactly the same as for the 7 r + 7 r - 7 
analysis. The remaining pion background after the 
p + p~'Y selection is much reduced, reaching 10 % in 
the p peak region. A comparison between data and 
MC is shown in Fig. [I] The same data sample as 
used in the main analysis is also used here, but we 
present a larger mass range than for the 7 t + 7 t - 7 
case. The efficiency corrections described in the 
previous section have been applied to MC on a 
track and photon candidate basis. The lower panel 
of Fig. |T] shows the relative discrepancy between 
data and MC. A good agreement over the full 
mass range at the level of (1.0 ± 0.3 ± 0.9)% and 
X 2 /ndf = 134/139 is found, where the uncertain¬ 
ties are statistical and systematic, respectively. A 
difference in the mass resolution due to detector 
effects between data and MC is visible around the 
narrow J/ip resonance. A fit in the mass range 600 
< TO/j/j, < 900 MeV/c 2 , which is the mass range 
studied in the main analysis, gives a relative dis¬ 
crepancy of (2.0 ± 1.7 ± 0.9)%; this is illustrated in 
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the inset of the upper panel of Fig. |TJ The theoret¬ 
ical uncertainty of the MC generator Phokhara is 
below 0.5% [16], while the systematic uncertainty of 
our measurement is 0.9%. The latter is dominated 
by the luminosity measurement, which is needed 
for the normalization of the data set. We consider 
the good agreement between the 7 QED pre¬ 
diction and data as a validation of the accuracy 
of our efficiency corrections. As a further cross 
check, we have applied the efficiency corrections 
also to a statistically independent /i + p “7 sample, 
resulting in a difference between data and MC of 
(0.7 ± 0.2)% over the full mass range, where the 
error is statistical only. 



Figure 1: Invariant n + n~ mass spectrum of data 
and 7 MC after using the ANN as muon se¬ 
lector and applying the efficiency corrections. The 
upper panel presents the absolute comparison of the 
number of events found in data and MC. The inset 
shows the zoom for invariant masses between 0.6 
and 0.9 GeV/c 2 . The MC sample is scaled to the 
luminosity of the data set. The lower plot shows the 
ratio of these two histograms. A linear fit is per¬ 
formed to quantify the data-MC difference, which 
gives a difference of (1.0 ± 0.3 ± 0.9)%. A differ¬ 
ence in the mass resolution between data and MC 
is visible around the narrow J/ip resonance. 


6. Extraction of er(e+e —> 7T+7T ) and |.F 2 | 

6.1. Methods 

We finally extract a „ = a(e + e~ —> 7r + 7r ~) ac¬ 
cording to two independent normalization schemes. 


In the first method, we obtain the bare cross sec¬ 
tion, i.e., the cross section corrected for vacuum 
polarization effects, according to the following for¬ 
mula: 


bare 

cr 7rir(7F S R) 


N, 


TVTVy 


C-e 


.7r7T'Y 

global 


(1 + *fsr) 

• H(s) ■ (5vac 


( 1 ) 


where N n7r ~ f is the number of signal events found 
in data after applying all selection requirements de¬ 
scribed above and an unfolding procedure to correct 
for the mass resolution, C the luminosity of the data 
set, and H the radiator function. The global effi¬ 
ciency determined based on the signal MC 

by dividing the measured number of events after all 
selection requirements AW“f sured by that of all gen¬ 
erated events Ag™(; rated . The true MC sample is 
used, with the full d 7 range, applying the efficiency 
corrections mentioned in Section 13.31 but without 
taking into account the detector resolution in the 
invariant mass m: 


^global (^) 


^meTsuredM 

derated M 


( 2 ) 


The efficiency is found to depend slightly on m nn 
and ranges from 2.8% to 3.0% from lowest to high¬ 
est m W7r . An unfolding procedure, which eliminates 
the effect of the detector resolution, is described in 
Sect. |6.2| and is applied before dividing by the global 
efficiency. The radiator function H is described in 


Sect. 6.4 As input for a„ the bare cross section is 


needed. It can be obtained by dividing the cross 
section by the vacuum polarization correction <5 vac , 


which is also described in Sect. |6.4| As pointed out 
in Ref. HU. in order to consider radiative effects 
in the dispersion integral for a M , an FSR correction 
has to be performed. The determination of the cor¬ 


rection factor (1 4 - <5Jg R ) is described in Sect. 6.3 


In the second method, we use a different nor¬ 
malization than in the first method and normalize 
N-k-kj to the measured number of /i + /x _ 7 events, 
Since C , H , and <5 vac cancel in this normal¬ 
ization, one finds the following formula: 


o. 


1\T 1 , 

bare _ . global . 1 ■ °FSR 

7T7t(7fSr) l\f 7T7T7 1 i ?7T7T 

^global 1 T °FSR 


where is the global efficiency of the dimuon 

selection, already described in Sect. 1 dpg R is the 
FSR correction factor to the final state, 

which can be obtained using the Phokhara event 
generator, cr^ re is the exact QED prediction of the 
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dimuon cross section, given by [2BJ Eq. (5.13)] 

bare M 3 ~ 

^ = ~ 3 ? - 2 - ’ (4) 

with the fine structure constant a, the cms en¬ 
ergy s' < s available for the creation of the final 
state, the muon velocity /3 M = y / T^4m 2 7s', and 
the muon mass m^. The contributions of radia¬ 
tor function, luminosity, and vacuum polarization 
to the systematic uncertainties of the bare cross 
section, cancel in the second method. The upper 
panel of Fig. [2] shows the comparison of the bare 
cross sections including FSR obtained with the first 
(black) and second method before unfolding (blue). 
The error bars are statistical only. They are much 
larger for the second method due to the limited 
p + p ~7 statistics in the mass range of interest. The 
lower panel shows the ratio of these cross sections. 
Again, a linear fit is performed to quantify the dif¬ 
ference, which is found to be (0.85 ± 1.68)% and 
X 2 /ndf = 50/60, where the error is statistical. Both 
methods agree within uncertainties. The first one 
is used in the analysis. Finally, the pion form factor 
as a function of s' can be calculated via 

l^lV) = 2oi n ^ eSSe V) , (5) 

7 rarp£(s') 

with the pion velocity /3 OT (s') = y/l — Am^/s', the 
charged pion mass m^, and the dressed cross sec¬ 
tion CT^ essed (s') = cr(e + e~ —> n + Tr~)(s') containing 
vacuum polarization, but corrected for FSR effects. 
The result is presented in Sect. [7] 



ifs 7 [GeV] 




** 


Figure 2: Comparison between the methods to ex¬ 
tract ov-n- explained in the text — using the lumi¬ 
nosity (black) and normalizing by (blue). The 
lower panel shows the ratio of these results together 
with a linear fit (blue line) to quantify their differ¬ 
ence. 


6.2. Unfolding 

In order to obtain the final result for 
one has to rectify the detector resolution effects, 
i.e., the mass spectrum needs to be unfolded. To 
this end, the Singular Value Decomposition (SVD) 
method 123 is used. It requires two input variables 

the response matrix and the regularization pa¬ 
rameter t. The SVD algorithm calculates an opera¬ 
tor which cancels the detector smearing by inverting 
the response matrix. We obtain the response ma¬ 
trix in the full mass range between threshold and 
3.0 GeV, using a signal MC sample. The matrix 
corresponds to the correlation of the reconstructed 
m nn spectrum, and the originally generated m W7r 
values. With the choice of a bin width of 5 MeV/c 2 , 
about 43% of events are found to be on the diagonal 
axis. 

To find the value of the regularization param¬ 
eter r, we compare two independent methods, as 
suggested in Ref. [27]. On the one hand, we per¬ 
form a MC simulation where r is optimized such 
that unfolded and true distributions have the best 
agreement. On the other hand, we process an al¬ 
gorithm, described in [27], exploiting the singular 
values of the response matrix. Both methods favor 
a similar regularization parameter of r = 72. 

To estimate the systematic uncertainties and to 
test the stability of the SVD method, we perform 
two cross checks. In both cases we use a 7r + 7r _ 7 MC 
sample which is independent of the one used to de¬ 
termine the response matrix. We modify and then 
unfold the spectra in both checks. In the first cross 
check, the reconstructed spectrum is smeared with 
an additional Gaussian error, which results in an 
about 20% larger detector smearing than expected 
from MC simulation. The resulting unfolded spec¬ 
trum reproduces the true one on the sub per mille 
level. In the second cross check, the mass of the 
p-resonance is varied systematically in the simula¬ 
tion in steps of 10 MeV/c 2 between 750 and 790 
MeV/c 2 . The response matrix is kept fixed and 
was determined with a p mass of 770 MeV/c 2 . In 
all cases, the masses of the p peak after unfold¬ 
ing are found to be close to the initially simulated 
masses. From the comparisons of these checks, we 
take the maximum deviation of 0.2% as systematic 
uncertainty. 

6.3. FSR correction 

The correction factor <5fsr is determined with 
the Phokhara generator in bins of m wn . Two dif- 




















ferent correction methods are used on the data to 
cross check whether it is applied correctly. 

(1) The whole FSR contribution of the 7r + 7r”7 
events is calculated with Phokhara, by dividing a 
true MC spectrum including FSR in NLO by the 
spectrum without any FSR contribution. The re¬ 
sulting distribution is used to correct data. As 
pointed out in Ref. m, for the dispersion inte¬ 
gral for a M , the FSR correction for the process 
e + e“ —► 7r + 7r~ needs then to be added again. We 
use the calculation by Schwinger assuming point¬ 
like pions: 


_dressed _ —dressed 

Cr 7T7I-(7) — a TTTT 


/ \ OL 

1 +v(s)~ 

7r 


( 6 ) 


where r](s) is the theoretical correction factor taken 
from [2S]. In the p-peak region it is between 0.4% 
and 0.9%. 

(2) A special version of the Phokhara generator 
is used [22|, which, in contrast to the standard ver¬ 
sion of the generator, distinguishes whether a pho¬ 
ton is emitted in the initial or the final state. In 
events in which photons have been radiated solely 
due to ISR, the momentum transfer of the virtual 
photon s 7 * is equal to the invariant mass of the two 
pions m 2 ^. However, if an FSR photon is emitted, 
the invariant mass is lowered due to this effect and 
hence < s 7 ». The effect can be removed by ap¬ 
plying an unfolding procedure, using again the SVD 
algorithm. Here, the response matrix is m 2 ^ vs. 
Sry *, obtained from a MC sample that includes FSR 
in NLO. The regularization parameter r is deter¬ 
mined as described in Sect. |6.2| After applying the 
corrections for the radiative 7r + 7r _ 7 process, which 
are of the order of 2%, one obtains the 7 t + 7t _ (7fsr) 
cross section directly. 

The difference between both methods is found 
to be (0.18 ± 0.13)%. Both methods are comple¬ 
mentary and agree with each other within errors. 
The difference is taken as systematic uncertainty. 
Finally, the correction obtained with method (1) is 
used in the analysis. 


6-4- Radiator function and vacuum polarisation 
correction 

The radiator function is implemented within 
the Phokhara event generator with NLO precision. 
Hence, a very precise description is available with 
a claimed uncertainty of 0.5% [I TS] . 

To obtain the bare cross section, vacuum po¬ 
larization effects <5 vac must be taken into account. 


To this aim, the dressed cross section, including 
the vacuum polarization effects, is adjusted for the 
running of the coupling constant a [30] . Bare and 
dressed cross sections are related as follows: 



The correction factors are taken from Ref. 311] ■ 

6.5. Summary of systematic uncertainties 

Systematic uncertainties are studied within 
the investigated m W7r range between 600 and 
900 MeV/c 2 . Sources are: 

(1) Efficiency corrections: Each individual uncer¬ 
tainty is studied in bins of m,„ with respect to 
three different sources. Firstly, the remaining back¬ 
ground contaminations in the data samples are esti¬ 
mated with the corresponding MC simulation men¬ 
tioned in Tab. 1. Their contribution is taken into 
account by multiplying the claimed uncertainties of 
the event generators and their fraction of the inves¬ 
tigated signal events. Secondly, we vary the selec¬ 
tion requirements {E/p, Xic- depth of a charged 
track in the MUC), which are used to select clean 
muon and pion samples for the efficiency studies, in 
a range of three times the resolution of the corre¬ 
sponding variable. The differences of the correction 
factors are calculated. Thirdly, the resolution of the 
correction factors, i.e., the bin sizes of momentum 
and 6 distributions, is varied by a factor two and 
the effects on the final correction factors are tested. 

(2) Pion-muon separation: Additional uncertainties 
of using the ANN method for pion-muon separa¬ 
tion are estimated by comparing the result from a 
different multivariate method, namely the Boosted 
Decision Tree (BDT) approach [23]. As a further 
cross check, the whole analysis is repeated without 
the use of a dedicated PID method. 

(3) Residual background is subtracted using simu¬ 
lated events. The uncertainty is determined to be 
0 . 1 %. 

(4) Angular acceptance: The knowledge of the an¬ 
gular acceptance of the tracks is studied by varying 
this requirement by more than three standard de¬ 
viations of the angular resolution and studying the 
corresponding difference in the selected number of 
events. A difference of 0.1% in the result can be 
observed. The procedure is repeated for all other 
selection criteria. Their contribution to the total 
systematic uncertainty is found to be negligible. 

(5) Unfolding: Uncertainties introduced by unfold¬ 
ing are smaller than 0.2%, as estimated by the two 
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cross checks mentioned in Sect. 16.21 

(6) FSR correction: The uncertainty due to the 
FSR correction is obtained by comparing two dif¬ 
ferent approaches as described in Sect. |6.3| The 
uncertainty is found to be 0.2%. 

(7) Vacuum Polarization: The uncertainty due to 
the vacuum polarization correction is conservatively 
estimated to be 0.2%. 

(8) Radiator Function: The Radiator Function 
extracted from the Phokhara generator is imple¬ 
mented with a precision of 0.5%. 

(9) Luminosity: The luminosity of the analyzed 
data set has been determined to a precision of 0.5%. 
All systematic uncertainties are summarized in 
Tab. [2j They are added in quadrature, and a 
total systematic uncertainty for cr bare (e + e~ —► 
7t + 7t - (7fsr)) of 0.9% is achieved, which is fully cor¬ 
related amongst all data points. 

Table 2: Summary of systematic uncertainties. 


Source 

Uncertainty 

(%) 

Photon efficiency correction 

0.2 

Pion tracking efficiency correction 

0.3 

Pion ANN efficiency correction 

0.2 

Pion e-PID efficiency correction 

0.2 

ANN 

negl. 

Angular acceptance 

0.1 

Background subtraction 

0.1 

Unfolding 

0.2 

FSR correction <5fsr 

0.2 

Vacuum polarization correction S vac 

0.2 

Radiator function 

0.5 

Luminosity C 

0.5 

Sum 

0.9 


7. Results 

The result for <r bare (e + e _ —► 7r + 7r _ (7 F sR)) as a 
function of ^fs = m„ is illustrated in Fig. [ 3 ] and 
given numerically in Tab. [4j The cross section is 
corrected for vacuum polarization effects and in¬ 
cludes final state radiation. Besides the dominant 
p(770) peak, the well-known structure of the p-ui in¬ 
terference is observed. The result for the pion form 
factor \F n \ 2 is shown in Fig.[4]and given numerically 
in Tab. [4] It includes vacuum polarization correc¬ 
tions, but, differently from the cross section shown 
in Fig. [3j final state radiation effects are excluded 
here. The red line in Fig. [4] illustrates a fit to data 


Table 3: Fit parameters and statistical errors of the 
Gounaris-Sakurai fit of the pion form factor. Also 
shown are the PDG 2014 values [55] . 


Parameter 

BESIII value 

PDG 2014 

Trip [MeV/c 2 ] 

776.0 ± 0.4 

775.26 ± 0.25 

T p [MeV] 

151.7 ± 0.7 

147.8 ± 0.9 

m u [MeV/c 2 ] 

782.2 ± 0.6 

782.65 ± 0.12 

[MeV] 

fixed to PDG 

8.49 ± 0.08 

M [ 10 - 3 ] 

1.7 ± 0.2 

- 

\<j>u\ [rad] 

0.04 ± 0.13 

- 



Figure 3: The measured bare e + e —> tt + tt (7fsr) 
cross section. Only the statistical errors are shown. 



Figure 4: The measured squared pion form factor 
| | 2 . Only statistical errors are shown. The red 
line represents the fit using the Gounaris-Sakurai 
parametrization. 

according to a parametrization proposed by 
Gounaris and Sakurai [32]. Here, exactly the same 
fit formula and fit procedure are applied as de¬ 
scribed in detail in Ref. m- Free parameters of 


10 

















the fit are the mass and width T of the p meson, 
the mass of the u> meson, and the phase of the Breit- 
Wigner function c u = . The width of the 

u> meson is fixed to the PDG value [33]. The re¬ 
sulting values are shown in Tab. [3j As can be seen, 
the resonance parameters are in agreement with the 
PDG values [23] within uncertainties, except for P p , 
which shows a 3.4<r deviation. Corresponding am¬ 
plitudes for the higher p states, p(1450), p(1700), 
and p(2150), as well as the masses and widths of 
those states were taken from Ref. uni, and the sys¬ 
tematic uncertainty in F p due to these assumptions 
has not been quantitatively evaluated. 

The Gounaris-Sakurai fit provides an excellent 
description of the BESIII data in the full mass range 
from 600 to 900 MeV/c 2 , resulting in x 2 /ndf = 
49.1/56. Figure [5] shows the difference between fit 
and data. Here the data points show the statistical 
uncertainties only, while the shaded error band of 
the fit shows the systematic uncertainty only. 

In order to compare the result with previous 
measurements, the relative difference of the BE¬ 
SIII fit and data from BaBar ITT)] . KLOE [6, |7J Bj, 
CMD2 HU!, and SND [3 is investigated. Such a 
comparison is complicated by the fact, that previ¬ 
ous measurements used different vacuum polariza¬ 
tion corrections. Therefore, we consistently used 
the vacuum polarization correction from Ref. m 
for all the comparisons discussed in this section. 
The KLOE 08, 10, 12, and BaBar spectra have, 
hence, been modified accordingly. The individual 
comparisons are illustrated in Figs. [5] and [ 6 ] Here, 
the shaded error band of the fit includes the sys¬ 
tematic error only, while the uncertainties of the 
data points include the sum of the statistical and 
systematic errors. We observe a very good agree¬ 
ment with the KLOE 08 and KLOE 12 data sets 
up to the mass range of the p-cu interference. In the 
same mass range the BaBar and KLOE 10 data sets 
show a systematic shift, however, the deviation is, 
not exceeding 1 to 2 standard deviations. At higher 
masses, the statistical error bars in the case of BE¬ 
SIII are relatively large, such that a comparison is 
not conclusive. There seem to be a good agreement 
with the BaBar data, while a large deviation with 
all three KLOE data sets is visible. There are indi¬ 
cations that the BESIII data and BESIII fit show 
some disagreement in the low mass and very high 
mass tails as well. We have also compared our re¬ 
sults in the p peak region with data from Novosi¬ 
birsk. At lower and higher masses, the statistical 
uncertainties of the Novosibirsk results are too large 


to draw definite conclusions. The spectra from SND 
and from the 2006 publication of CMD-2 are found 
to be in very good agreement with BESIII in the p 
peak region, while the 2004 result of CMD-2 shows 
a systematic deviation of a few percent. 

We also compute the contribution of 
our BESIII cross section measurement 
crbare^ e + e — n + w - (y FSR )) to the hadronic con¬ 
tribution of (g — 2) p , 

(0.9GeV) 2 

a p 7r ’ LO (0.6 - 0.9 GeV) = ^ J ds'K^')^ , 

(0.6GeV) 2 

(8) 

where K(s') is the kernel function [IT] Eq. (5)]. 
As summarized in Fig. [7] the BESIII result, 
a ^,L°(6oo-900 MeV) = (368.2 ± 2.5 sta t ± 3.3 sys ) • 
10 ” 10 , is found to be in good agreement with all 
three KLOE values. A difference of about 1.7er with 
respect to the BaBar result is observed. 

8. Conclusion 

A new measurement of the cross section 
ffbare^g+g— 7r+7r _ ( 7 FSR )) tias been performed 
with an accuracy of 0.9% in the dominant p(770) 
mass region between 600 and 900 MeV/c 2 , using 
the ISR method at BESIII. The energy dependence 
of the cross section appears compatible with cor¬ 
responding measurements from KLOE and BaBar 
within approximately one standard deviation. The 
two-pion contribution to the hadronic vacuum po¬ 
larization contribution to (g — 2) p has been deter¬ 
mined from the BESIII data to be a p 7r,LO (600 — 
900 MeV) = (368.2 ± 2.5 sta t ± 3.3 sys ) • 10" 10 . By 
averaging the KLOE, BaBar, and BESIII values of 
a 7T7i-,L° anc j assum i n g that the five data sets are in¬ 
dependent, a deviation of more than 3cr between the 
SM prediction of (g — 2) p and its direct measure¬ 
ment is confirmed. For the low mass region < 600 
MeV/c 2 and the high mass region > 900 MeV/c 2 , 
the BaBar data was used in this calculation. 
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Figure 6: Relative difference of the form factor 
squared from KLOE 0 [7J [8] and the BESIII 
fit. Statistical and systematic uncertainties are in¬ 
cluded in the data points. The width of the BESIII 
band shows the systematic uncertainty only. 
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Table 4: Results of the BESIII measurement of the cross section er b + r ®_ j 7fsr ^ = cr bare (e + e —> tt + tt (tfsr)) 

and the squared pion form factor | | 2 . The errors are statistical only. The value of Vs' represents the bin 

center. The 0.9% systematic uncertainty is fully correlated between any two bins. 


%/s 7 [MeV] 

a ^-(7 F sa) t nb l 

l^l 2 

Vs' [MeV] 

^-(-vfsk) t nb ] 

l^l 2 

602.5 

288.3 ± 15.2 

6.9 ± 0.4 

752.5 

1276.1 ± 29.8 

41.8 ± 1.0 

607.5 

306.6 ± 15.5 

7.4 ± 0.4 

757.5 

1315.9 ± 31.3 

43.6 ± 1.0 

612.5 

332.8 ± 16.3 

8.2 ± 0.4 

762.5 

1339.3 ± 30.9 

44.8 ± 1.0 

617.5 

352.5 ± 16.3 

8.7 ± 0.4 

767.5 

1331.9 ±: 30.8 

45.0 ± 1.0 

622.5 

367.7 ± 16.6 

9.2 ± 0.4 

772.5 

1327.0 ± 30.6 

45.2 ± 1.0 

627.5 

390.1 ± 17.7 

9.8 ± 0.4 

777.5 

1272.7 ± 29.2 

43.7 ± 1.0 

632.5 

408.0 ± 18.0 

10.4 ± 0.5 

782.5 

1031.5 ± 26.7 

37.1 ± 0.9 

637.5 

426.6 ± 18.1 

11.0 ± 0.5 

787.5 

810.7 ± 24.2 

30.3 ± 0.8 

642.5 

453.5 ± 19.0 

11.8 ± 0.5 

792.5 

819.7 ± 23.8 

30.6 ± 0.8 

647.5 

477.7 ± 18.5 

12.5 ± 0.5 

797.5 

803.1 ± 23.3 

30.1 ± 0.8 

652.5 

497.4 ± 19.5 

13.2 ± 0.5 

802.5 

732.4 ± 22.1 

27.7 ± 0.8 

657.5 

509.2 ± 19.4 

13.6 ± 0.5 

807.5 

679.9 ± 20.6 

25.9 ± 0.7 

662.5 

543.4 ± 19.9 

14.7 ± 0.5 

812.5 

663.6 ± 21.0 

25.5 ± 0.8 

667.5 

585.0 ± 20.5 

16.0 ± 0.6 

817.5 

622.2 ± 19.9 

24.1 ± 0.7 

672.5 

642.7 ± 22.2 

17.7 ± 0.6 

822.5 

585.0 ± 19.5 

22.9 ± 0.7 

677.5 

640.5 ± 21.0 

17.8 ± 0.6 

827.5 

540.8 ± 18.1 

21.4 ± 0.7 

682.5 

668.0 ± 21.9 

18.8 ± 0.6 

832.5 

496.4 ± 17.7 

19.8 ± 0.7 

687.5 

724.4 ± 22.9 

20.6 ± 0.6 

837.5 

450.4 ± 16.8 

18.1 ± 0.6 

692.5 

783.5 ± 23.2 

22.5 ± 0.7 

842.5 

404.7 ± 15.2 

16.4 ± 0.6 

697.5 

858.6 ± 25.3 

24.9 ± 0.7 

847.5 

391.3 ± 15.4 

16.0 ± 0.6 

702.5 

893.8 ± 25.4 

26.2 ± 0.7 

852.5 

364.0 ± 15.0 

15.0 ± 0.6 

707.5 

897.8 ± 25.0 

26.6 ± 0.7 

857.5 

339.6 ± 14.0 

14.2 ± 0.6 

712.5 

978.6 ± 26.6 

29.3 ± 0.8 

862.5 

310.0 ± 13.7 

13.0 ± 0.6 

717.5 

1059.1 ± 27.9 

32.0 ± 0.8 

867.5 

283.8 ± 13.0 

12.1 ± 0.5 

722.5 

1086.0 ± 28.3 

33.2 ± 0.9 

872.5 

256.5 ± 12.4 

11.0 ± 0.5 

727.5 

1088.4 ± 27.7 

33.6 ± 0.9 

877.5 

237.3 ± 11.4 

10.3 ± 0.5 

732.5 

1158.8 ± 29.2 

36.2 ± 0.9 

882.5 

229.7 ± 11.6 

10.0 ± 0.5 

737.5 

1206.5 ± 29.6 

38.2 ± 0.9 

887.5 

224.0 ± 11.6 

9.9 ± 0.5 

742.5 

1229.9 ± 29.0 

39.3 ± 0.9 

892.5 

196.1 ± 10.5 

8.7 ± 0.4 

747.5 

1263.3 ± 30.3 

40.9 ± 1.0 

897.5 

175.9 ± 9.7 

7.9 ± 0.4 
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